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Abstract 

We analyze stability of a thin inextensible elastic rod which has non-vanishing 
spontaneous generalized torsions in its stress-free state. Two classical problems are 
studied, both involving spontaneously twisted rods: a rectilinear rod compressed by 
axial forces, and a planar circular ring subjected to uniform radial pressure on its outer 
perimeter. It is demonstrated that while spontaneous twist stabilizes a rectilinear rod 
against buckling, its presence has an opposite effect on a closed ring. 

The history of the buckling instability in a thin elastic rod under compression goes back to 
the works of Euler and his contemporaries [|l| , who laid down the foundations of the general 
theory of elastic stability 0. Although this field has long been the domain of engineers and 
applied mathematicians, there has recently been a renaissance of interest in problems related 
to the stability of thin filaments in the theoretical physics community 0-[16| prompted by 



experimental advances in the art of mechanical manipulation of single DNA molecules and 



rodlike protein assemblies |T^-[|l^. Because double stranded DNA molecules are helices 



that contain strongly curved regions |ED[, their modeling as elastic rods has to take into 



account their spontaneous curvature and twist. While general considerations concerning 
such "naturally" curved rods can already be found in the work of Kirchhoff little is 
known about the effect of spontaneous curvature and twist on the stability with respect to 
buckling. The present study deals with the effect of spontaneous twist on the static stability 
of a thin inextensible elastic rod. Nonlinear equilibrium equations are derived for a rod under 
the action of external forces and moments. Two classical problems are analyzed which have 
well-known solutions in the absence of spontaneous twist a rectilinear rod compressed 
by forces applied to its ends and a ring compressed by radially inwards directed pressure 
uniformly distributed along its perimeter. It is demonstrated that while spontaneous twist 
stabilizes the straight rod against buckling, it has a destabilizing effect on the ring. In both 
cases buckling takes place through a three-dimensional instability. 



1 Equilibrium equations 

Consider an elastic rod of length L and cross-section S which has two axes of symmetry. We 
assume that a where a is the largest dimension associated with the cross-section. Denote 



1 



by s the arc-length which parameterizes the centerhne of the rod in the stress-free reference 
configuration. The radius vector Ro(s) describes the spatial position of the centerline and 
to = dHo/ds the unit tangent vector in the reference state. The unit vectors uo(s) and vo(s) 
are directed along the axes of symmetry of the cross-section. We would like to emphasize 
that Uq and vq are associated with physical attributes of rods of small but non-vanishing 
thickness and, in general, do not coincide with the normal and the bi-normal vectors which 
appear in the description of three-dimensional configurations of infinitely thin geometrical 
lines. Differentiation of the vectors tg, uq and vq with respect to s results in 

dto duo dvo 

ds ds ds 

where uJko{s) are generalized intrinsic torsions. The position of an arbitrary point of the rod 
is given by the Lagrangian coordinates {^fc}, where ^1,^2 are Cartesian coordinates in the 
cross-sectional plane with unit vectors uq and vq, and ,^3 = s, 

ro(^i, 6, s) = Ro(s) + ^iuo(s) + 6vo(s). 

We assume that the centerline is inextensible, which means that s remains the arc-length in 
the deformed state as well. The position of the centerline in the deformed configuration is 
determined by its radius vector R(s). The position of an arbitrary point of the cross-section 
with the arc-length s is 

r(6, 6, s) = R{s) + ^^u{s) + 6v(s), (2) 

where u and v are unit vectors directed along the axes of symmetry of the cross-section, 
and t = dH/ds is the unit tangent vector. The functional form of Eq. is based on the 
assumption that after deformation any cross-section remains planar and perpendicular to 
the centerline of the rod. By analogy with Eq. (|I]), we write 

dt du dv 

— = IViU + UJ2V, — = -UJit + CUsV, — = -Co'2t - UJ3U, (3) 

ds ds ds 

where u!k{s) are generalized torsions. 



Using the methods of Ref. ||2T[, we calculate the displacement gradient in the deformed 
state dro/dr and the strain tensor. The mechanical energy is constructed as a quadratic 
form in this strain tensor. Expressions for the components of the stress tensor are derived 
by differentiating this energy with respect to the corresponding components of the strain 
tensor. The constitutive equations relate the components of the internal moment M = 
M^u + MyV + Mft to deviations of the generalized torsions from their spontaneous values, 

Mu = -A2 (W2 - CJ20) , M^ = Ai {ui - cjio) , Mt = A3 (W3 - W30) , (4) 

where Ai = Eili, A^ = E1I2, A3 = E2{Ii + 12)1 h and I2 are moments of inertia with 
respect to the two symmetry axes of the cross-section, and Ei and E2 are elastic moduli. 
Denote the internal force by F = EuU + F„v + Ftt, the external distributed force per unit 
length by q = g„u -|- g„v + qtt, and the external distributed moment per unit length by 
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m = m^u + m^jV + m^t. Upon some algebra, we eliminate the components of the internal 
force in the cross-sectional plane and write down the conditions of mechanical equilibrium, 

dFt dM^ dMu , , , 

+ ui— UJ2—, h uj^[uJiMu + uj2M.^) + uJim^ - uj2mu + = 0, 

ds ds ds 

- UiMu - uj2M^ + mt = 0, 

ds 

d . dMy dMu , ^ , , N 

— (— h uj^Mu + uj2Mt + m^) + UJ2,[— uj^M-u + UiMt + m„) 

ds ds ds 

-UiFt -qu = 0, 

d , dMu dMu 

— (— uj^Mu + iOiMt + niu) - t^3(— 1 \- ^zMu + uj2Mt + m„) 

ds ds ds 

+uj2Ft + = 0. (5) 

Given boundary conditions, the four nonlinear differential equations (|^) together with the 
constitutive relation (H), uniquely determine the four unknown functions, F^, oji, uj2 and u^. 
The generalized torsions ujk are substituted into Eqs. (|^) which are resolved to determine 
the spatial configuration of the centerline and the rotation of the cross-section (twist) of the 
deformed rod. The tangential force Ft plays the role of a Lagrange multiplier which enforces 
the inextensibility condition. 



2 Stability of a rectilinear rod 

We begin with the analysis of stability of a spontaneously twisted rectilinear rod. In the 
stress-free reference state, the rod's centerline is straight, and its cross-sections rotate around 
the centerline by some angle Vt = Q{s). We set Ai < A2, which means that the vector uq 
coincides with the smaller of the two principal axes of inertia. The vectors uq, vq and to 
are given by uq = cosf2ei + sinne2, vq = — sinfiei + cosf2e2, and tg = 63, where are 
unit vectors of some Cartesian coordinate frame associated with the stress-free rod. These 
equalities and Eq. imply that 

dn 

uio = 0, UJ20 = 0, UJ30 = —. 

ds 

We will assume in the following that u^q is independent of s. 

When some compressive force P (smaller than the critical force for buckling) is applied 
to the ends of an inextensible rod, the rod does not deform, the moments Mu, M^ and 
Mt vanish, and only a longitudinal internal force —P is generated. We now assume that, 
in addition to the load P, small perturbations of external distributed forces and external 
distributed moments are present. These perturbations cause small displacements of the rod 
that consist of small displacements of the centerline, AR = Xiei + ^262, and small rotations 
of the cross-sections around the centerhne a. Under this perturbation, the vectors uq, vq and 
to are replaced by u, v, t, which are substituted into Eq. to obtain expressions for the 
perturbed generalized torsions Uk (linear in the perturbations xi, X2 and a). Substitution of 
these expressions into Eq. (^ yields the perturbations of internal moments, 

. /(i^xi . ^ d'^X2 ^\ , /(i^Xi ^ d'^X2 . ^\ .da 

Mu = A2[-T^smn- —^cosn), Mu = Ai[—-^cosn + —^smn), Mt = A^ — . 
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Substituting these expressions into Eqs. we find that the equations for the perturbation 
of the longitudinal force and for the rotation a decouple from the equations for the curvatures 
io\ and ijj2i which govern bending of the rod and determine it stability with respect to 
buckling. 

When the external load P reaches its critical value, the rod becomes unstable against 
arbitrarily small external perturbations of moments and forces and, therefore, buckling can 
be interpreted as the point at which a non- vanishing solution to Eqs. (j^), in which these 
forces and moments are set to zero, first appears. This yields 

'~d^ + - ^1^30)^l - "^^2^30^ = 0, 
A2^ + (P - A2iulo)^2 + 2^1^30^ = 0. (6) 

Equations (H) are supplemented by the geometric relations which follow from Eqs. (|^): 

d'^Xi d'^X2 „ „ 

= LJi cos il — LJ2 sm i 2, = uji sm il + uj2 cos i I. 

ds^ ds^ 

The resulting linear equations are solved subject to the boundary conditions corresponding 
to hinged ends: 

xi(0) = x,{L) = X2{0) = X2{L) = 0, ^(0) = ^(L) = ^(0) = ^(L) = 0. 

ds^ as^ as^ cts^ 

When the spontaneous twist u^q vanishes, Eqs. turn into conventional equations for 
stability of an elastic beam 0. In the general case, Eqs. (H) may be resolved with respect 
to Ui and uj2- Each of these functions obeys a linear fourth-order differential equation which 
contains only even derivatives with respect to s and can, therefore, be written as a linear 
combination of sm{j3ks) and cos(/5fcs), where the "wavenumbers" jSi and (32 are the positive 
roots of a characteristic equation that guarantees the existence of non-vanishing solutions 
uji{s) and 002(3) to Eqs. (||). Using standard methods, we derive a transcendental equation 
which relates the dimensionless critical force p* = Per /Per, where = vr^Ai/L^ is the Euler 
critical force for an untwisted rod [|l], to the total spontaneous twist angle Q{L) = uj^qL. 
The numerical solutions of these equations for various values of the cross-section asymmetry 
u = I2/I1 are plotted in Figure 1. Note that at the point of the instability, all generalized 
torsions Uk do not vanish, which implies that, unlike the planar Euler instability, buckling 
of a spontaneously twisted rod is intrinsically three-dimensional. 

For a circularly symmetric cross-section (z/ = 1), the Euler result holds for arbitrary 
values of the spontaneous twist, i.e., p* is independent of fl. For any u > 1, the critical 
force increases initially with Q{L) and reaches a maximum at Q{L) = 27t. Upon further 
increase in the total angle of twist, the critical force oscillates (with additional maxima at 
integer multiples of vr) and eventually tends to some limiting value which depends on u. 
To understand the origin of this non-monotonic behavior, we investigated the behavior of 
critical wavenumbers f3k- Numerical calculations show that the difference between the critical 
wavenumbers at buckling, (32 — Pi, oscillates as a function of ^1{L), with maxima occurring at 
integer multiples of vr, and tends asymptotically to 27r. In this asymptotic limit the critical 
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wavenumbers approach the "resonance" conditions l3kL — >• ui^qL^t: {k = 1,2). We conclude 
that the oscillatory dependence of the critical force on the twist angle arises as a result of 
the interplay between the two modes. Whenever the total angle of twist coincides with an 
integer multiple of vr, the corresponding axes of symmetry of the cross-sections at both ends 
of the rod (at s = and s = L) become aligned, just as in the twist-free case. 



3 Stability of a ring 

We now turn to study the stability of a circular ring whose centerline in the stress-free state 
is a planar circle with radius Rq- We introduce cylindrical coordinates {r, 0, z} with unit 
vectors e^, e^, e^, where the axis z is perpendicular to the plane of the ring. The initial 
position of the centerline in the stress-free reference state is determined by the radius vector 
Ro = RoSr. We assume that in its stress free reference state the rod is characterized by a 
twist angle Q = Q{s) between the vectors — and uq. The unit vectors to, Uq and vq can 
be expressed in terms of e^, and as 

to = e,^, Uo = — cos flSr + sin Q e^, vo = sin Q e^. + cos Q e^. 

Substitution of these expressions into Eqs. (|lD results in 

1 1 . dn 

ujio = —cosil, u;20 = -^smi2, uj3o = —. 
rCo iXQ ds 

We study the stability of the twisted ring under a radially inward-directed force per unit 
length q = — ge^ that acts on the outer perimeter of the ring. As long as the force is smaller 
than some critical value, it is balanced by a tangential internal force Ft = —qRo and no 
displacement arises. We consider a small deviation of the centerline from its unperturbed 
position, Ro ^ R = Ro + yr^r + y<j,G^ + Vz^z, and a small perturbation of the rotation 
angle about the centerline, Q —>■ Q + a. We express the perturbed unit vectors t, u and 
V in terms of the small deviations y,^, i/z and a {nr is eliminated using the inextensibility 
condition |t| = 1) and use Eq. to relate these deviations to the perturbed spontaneous 
torsions Uk- The subsequent analysis is qualitatively similar (but much more cumbersome) 
to that leading to Eqs. @. We derive three coupled differential equations with variable 
coefficients for the functions y(j,{s), yz{s) and a{s). An analytical solution of these equations 
can be obtained only in the special case of equal moments of inertia, Ji = I2 (e.g., for a 
circular cross-section). In this case we find that buckling takes place when the critical force 
per unit length reaches a critical value 

_ SE,h 

As expected for a circular rod, the critical force does not depend on the spontaneous twist 
Q. In this case, the onset of instability with respect to in-plane buckling of a circular rod 
coincides with that for out-of-plane buckling. 

In the general case, Ii < I2, we have to resort to approximate methods. For simplicity, we 
set Q{s) = Ns/Rq, where is a positive integer. A reasonable estimate which is expected 
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to capture the qualitative aspects of the effect of spontaneous twist on the critical force, can 
be obtained using the Galerkin method 0. The perturbations y^, Uz and a are presented in 
the form X„sin(ns/-Ro) with a fixed n. We substitute these expressions into the equilibrium 
equations, multiply both sides by sin(ns/-Ro) and require that the difference between the 
left and the right hand sides of the resulting equalities vanishes on the average (i.e., when 
integrated over the ring contour, from s = to 2ttRq). Upon some algebra, we arrive at the 
following expression for the critical compressive force per unit length 

Note that for N"^ > 3(/i + l2)/{h — h) — ^ the ring becomes unstable with respect to 
buckling, even in the absence of external forces. Because the equations for the perturbations 
are coupled, the instability has a three dimensional character and buckling is not confined 
to the plane of the ring. 

It is of interest to compare these conclusions with the results of Ref. where the 
stability of the stress-free configuration of a circular ring was analyzed and it was found 
that, in general, there were no non-trivial solutions to the equilibrium equations with periodic 
boundary conditions. The apparent discrepancy stems from the fact that the problems under 
consideration are rather different: while we examine the stability of a ring with spontaneous 
twist under compression in the plane of the ring, Ref. analyzes the stability of a ring 
without spontaneous twist under the action of an external torque applied to a fixed cross- 
section. 

4 Conclusions 

Buckling instabilities have been studied for thin rods with spontaneous twist and non-circular 
cross-sections. Two cases are analyzed in detail: (i) a rectilinear rod compressed by forces 
applied to its ends, and (ii) a circular ring compressed by uniformly distributed pressure 
acting in its plane. It is demonstrated that for a rectilinear rod, an increase in the spon- 
taneous twist angle leads to the growth of the critical force and stabilizes the rod against 
buckling. The effect increases with the ratio of principal moments of inertia, u = /2//1 > 1. 
For a circular ring, an increase in the spontaneous twist leads to a monotonic decrease in 
the critical force. The rate of the decrease grows with u. For sufficiently large values of the 
spontaneous twist, the ring may become unstable to buckling even in the absence of external 
forces. 

Acknowledgement 

We would like to thank S. Panyukov for helpful comments. AD gratefully acknowledges 
financial support by the Israeli Ministry of Science through grant 1202-1-98. YR's work is 
supported by a grant from the Israel Science Foundation. 

References 



6 



[1] Love, A.E.H. A Treatise on the Mathematical Theory of Elasticity (Dover, New York, 
1944). 

[2] Timoshenko, S.P. and Gere, J.M. Theory of Elastic Stability (McGraw-Hill, New York, 
1963). 

[3] Guitter, E. and Leibler, S. Europhys. Lett. 17 (1992) 643. 

[4] Marko, J.F. and Siggia, E.D. Macromolecules 27 (1994) 981. 

[5] MacKintosh, F.C., Kas, J. and Janmey, P.A. Phys. Rev. Lett. 75 (1995) 4425. 

[6] Goriely, A. and Tabor, M. Physica D 105 (1997) 20, 45. 

[7] Kamien, T.D., Lubensky, T.C., Nelson, P. and O'Hern, C.S. Europhys. Lett. 38 (1997) 
237. 

[8] Bouchiat, C. and M. Mczard, M. Phys. Rev. Lett. 80 (1998) 1556. 

[9] Goldstein, R.E., Powers, T.R. and Wiggins, C.H. Phys. Rev. Lett. 80 (1998) 5232. 

[10] Liverpool, T.B., Golestanian, R. and Kremer, K. Phys. Rev. Lett. 80 (1998) 405. 

[11] Marko, J.F. Phys. Rev. E 57 (1998) 2134. 

[12] Moroz, J.D. and Nelson, P. Macromolecules 31 (1998) 6333. 

[13] Balaeff, A., Mahadevan, L. and Schulten, K. Phys. Rev. Lett. 83 (1999) 4900. 

[14] Haijun, Z. and Zhong-can, O.-Y. J. Chem. Phys. 110 (1999) 1247. 

[15] Tobias, I., Swigon, D. and Coleman, B.D. Phys. Rev. E 61 (2000) 747, 759. 

[16] Garrivier, D. and Fourcade, B. Europhys. Lett. 49 (2000) 390. 

[17] Smith, S.B., Cui, Y. and Bustamante, C. Science 271 (1996) 795. 

[18] Strick, T.R., AUemand, J.-F., Bensimon, D. and Croquette, V. Science271 (1996) 1835. 

[19] Fygenson, D.K., Elbaum, M., Shraiman, B. and Libchaber, A. Phys. Rev. £'55 (1997) 
850. 

[20] Bednar, J., Furrer, P., Katritch, V., Statiak, A., Dubochet, J. and Stasiak, A. J. Mol. 
Biol. 254 (1995) 579. 

[21] Drozdov, A.D. and Rabin, Y. Phys. Rev. E (submitted). 



7 



TT 2n 3% 

1.25 1 1 1 r 




0.0 n{L) 10.0 



Figure 1: The critical load p* versus the total twist angle Q{L). Curve 1: v — 1.1; curve 2: 
u — 1.3; curve 3: u — 1.5 
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